Subgap tunneling via quantum-interference effect: insulators and charge density waves 
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A quantum interference effect is discussed for subgap tunneling over a distance comparable to the 
coherence length, which is a consequence of "advanced-advanced" and "retarded-retarded" trans- 
mission modes [Altland and Zirnbauer, Phys. Rev. B 55, 1142 (1997)]. Effects typical of disorder 
are obtained from the interplay between multichannel averaging and higher order processes in the 
tunnel amplitudes. Quantum interference effects similar to those occurring in normal tunnel junc- 
tions explain magnetoresistance oscillations of a CDWpierced by nanoholes [Latyshev et al., Phys. 
Rev. Lett. 78, 919 (1997)], having periodicity h/2e as a function of the flux enclosed in the 
nanohole. Subgap tunneling is coupled to the sliding motion by charge accumulation in the inter- 
rupted chains. The effect is within the same trend as random matrix theory for normal metal-CDW 
hybrids [Visscher et al, Phys. Rev. B 62, 6873 (2000)]. We suggest that the experiment by Laty- 
shev et al. probes weak localization-like properties of evanescent quasiparticles, not an interference 
effect related to the quantum mechanical ground state. 

PACS numbers: 73.20.Fz,73.23.-b,71.45.Lr,74.78.Na 



I. INTRODUCTION 

Normal electron tunneling 1 through a barrier is real- 
ized in two-terminal devices at the tip of a scanning tun- 
neling microscope, or by extended interfaces for planar 
tunnel junctions. Diffusive conductors can be described 
by arrays of tunnel junctions 2 ^ and it is desirable to de- 
velop a thorough understanding of all aspects of weak 
localization already at the level of a single tunnel junc- 
tion. Starting from a model system of normal electron 
tunneling through a band insulator (with hypothesis on 
the dimensionality and on the extrema of the disper- 
sion relation), we show that weak localization- like subgap 
tunneling 4 -! 5 - is the clue to an experiment by Latyshev et 
al& on h/2e oscillations of the magnetoresistance related 
to CDW motion around a nanohole of size comparable 
to the CDW coherence length. 

CDWs realize a well known phase of quasi-one dimen- 
sional (quasi-lD) conductors 7 -, with density modulations 
along the direction of the chains in the ground state, with 
a gap and a collective sliding motion above the depining 
threshold. CDWs can be nano-fabricated, as shown by 
various experiments in the last decade ^ 10 ' 11 ! 12 ' 13 ! 14 . 

A film of the CDW compound NbSe3 pierced by the 
nanoholes formed by columnar defects was obtained from 
ion irradiation by Latyshev et al.—, and the experimental 
results were reproduced after publication. The diame- 
ter D ~ 10 nm of the nanohole is comparable to the 
ballistic CDW coherence length £o = Hvp/A, with vf 
the Fermi velocity and A the Peierls gap of the CDW 6 . 
We reach an agreement with the following experimental 
observations 6 - for the (un) irradiated sample (not) con- 
taining nanoholes: 1. Absence of magnetoresistance os- 
cillations without nanoholes; 2. Absence of magnetoresis- 
tance oscillations with nanoholes but without sliding mo- 
tion; 3. h/2e oscillations of the resistance as a function of 
magnetic flux with nanoholes and with sliding motion; 4. 
Positive magnetoresistance at low fields with nanoholes; 
5. Oscillations measured by Latyshev et al& at temper- 



ature as high as ~ 52K in the presence of nanoholes. 

Previous approaches to related phenomena in CDWs 
were based on weak localization in normal metal-CDW 
hybrids in the framework of random matrix theory-i 5 -, 
Aharonov-Bohm oscillations 1 - 6 -, soliton tunnelin g 16 i 17 
and permanent current s 18 ! 19 . The effect that we consider 
is not directly related to non linearities of the CDW 
phase Hamiltonian 20 . 21 ! 22 ! 23 ! 24 . 25 ! 26 ! 27 ! 28 ! 29 ! 30 ! 31 ! 32 ^ 3 .. 
We reach consistency with Visscher et al. 15 finding 
on the basis of random matrix theory an unexpected 
strong effect of weak localization for CDWs con- 
nected to a disordered normal electrode. The effects 
discussed below rely on a geometrical parameter 
being comparable to the coherence length, namely 
D ~ £o for a nanohole of diameter D in a CDW 
film. By analogy, transport in superconducting hybrids 
having a dimension comparable to the superconduct- 
ing coherence length has aroused considerable interest 

rcccm -l „34, 35. 36, 37. 38, 39, 40,41,42, 43,44,45,46, 47,48. 49, 50. 51. 52, 53. 54, 55. 56 

in connection with the realization of a source of entangled 
pairs of electrons. 

Non standard localization effects were already intro- 
duced in superconducting hybrid structures by Altland 
and Zirnbauer for an Andreev quantum dot 5 - 7 -, who 
state at the end of their abstract that in normal metal- 
superconductor hybrids, "every Cooperon and diffusion 
mode in the advanced-retarded channel entails a corre- 
sponding mode in the advanced- advanced (or retarded- 
retarded) channel", which constitutes the technical basis 
of part of our discussion where "diffusion modes" become 
"transmission modes" . 

The structure of the article is as follows. Preliminaries 
are presented in Sec. [II] Sec. IIIII is focused on tunneling 
across band insulators. Sec. IIVI presents our results on 
charge density waves in connection with Sec. IIIII and with 
the experiments by Latyshev et aL— . Concluding remarks 
are presented in Sec. El The discussion in the main body 
of the article is complemented by two Appendices. 
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FIG. 1: (Color online.) Schematic representation of the BCS-like dispersion relation of the band insulator with a gap A. All 
states are populated by genuine electrons and holes, resulting in a band insulator, not in a superconductor. Electrons with 
wave-vectors ~ ±/cf are labeled by R and L respectively. A R electron can have a positive or negative group velocity, as shown 
on the figure. 



II. PRELIMINARIES: HAMILTONIANS, 
GREEN'S FUNCTIONS, WEAK LOCALIZATION 
LOOPS 



We start with necessary preliminaries on the micro- 
scopic description. Insulators (Sec. Ill A)) and charge den- 
sity waves (Sec. Ill B|) are described by the same formal- 
ism. Weak localization-like subgap tunneling is intro- 
duced in Sec. Ill CI 



A. Description of a band insulator 

We consider for the band insulator the same dispersion 
relation as for the BCS model (see Fig. [I}, but assume 
excitation branches supporting either electrons or holes, 
without condensate and without explicit CDW or super- 
conducting symmetry breaking. 

The advanced Green's function connecting two sites a 
and (3 separated by a distance R a ,p is given by 



sin(fcF|i? Q ,/3|) X 



(1) 



-Hu> 



: exp (-\R a ,p\/£(uj)), 



where pjq is the normal density of states, kp the Fermi 
wave-vector corresponding to th e extrcma of the disper- 
sion relation, and £(lu) = Hvf/ \/\A\ 2 — (hu>) 2 the coher- 
ence length at energy ftw, with vp the Fermi velocity. 

In addition, the insulator will be supposed to be quasi- 
1D in the forthcoming Sec. lIII Cl Then, we separate right 
(label R) from left (label L) branches (underlines are 
used to avoid confusion with the "A" and "R" labels 
for advanced and retarded Green's functions) having a 
wave- vector ~ ±kp. An electron on the R branch can 
have a positive or negative group velocity (see Fig. [TJ. 
The corresponding 2x2 matrix Green's function between 



times t and t' and positions a and {3 are given by 5? 

9aAt,t') = -i9(t-t')x (2) 



it'),Cp,L{t) 



After Fourier transform, Eq. © reduces to 

- A , s —Tijjt) 



exp(-|i?«,/3|/0 (3) 



VI A I 2 -(M 2 

exp {ikpRa^p) 

exp {—ikpR a ^p 



for propagation at energy hw between two points sep- 
arated by R a .j3 along a given chain. Eq. ^ in 
ID is compatible with Eq. (Q]) in 3D because the 
sin {kp\R a ^\)/kp\R a ^\ factor in 3D is replaced by 
cos (kpRa^) in ID. 

A coupling to the vector potential due to a magnetic 
field is obtained by replacing g£ Juj) in Eq. ([2]) by 



3a,/3(w)exp i- 



(4) 



where A • dr is the circulation of the vector potential 
on a path connecting a to f3, and 4>o = h/e is the flux 
quantum. 



B. Description of a charge density wave 

1. Peierls Hamiltonian 

Charge density waves are described on the basis of 
the electronic part of the Peierls Hamiltonian of spinless 
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(d) Elastic cotunneling in a charge density wave 

FIG. 2: (Color online.) Schematic representation in real space 
of a localization loop inserted in diagrams for the transmission 
coefficient in a normal metal (a); and in a gapped system 
below the gap (b). Similar diagrams are obtained from the 
tunnel terms at the interfaces^ to higher orders, (c) shows 
a weak localization-like loop for a quasi-lD charge density 
wave chains pierced by a nanohole (see Sec. IIV|I . (d) shows 
a process now winding around the hole, and therefore not 
being modulated by a magnetic flux in the nanohole. The 
labels "A" and "R" stand for advanced and retarded Green's 
functions. 



fermions on a ID chain: 



C x + C x Ca;4. O0 J 



(5) 

where the average hopping amplitude T has a modulation 
So (do is the lattice parameter). Changing variable from 
the spatial coordinate x along the chain to the wave- 
vector k, we find 



H = - 2T cos (*oo)Cfc c * 



(G) 



E Ac t-2k F C k - E A *4+2k F C k, 



k k 

with |A| = 5q the Peierls gap. 



B = f (u) exp (ik F (x a +xp)) (9) 

C = f (u)exp(-ik F (x a +xp)) (10) 

D = g (u) exp (~ik F (x a - xp)), (11) 

where the points a and (3 are at coordinates x a and xp 
along the chains, and where 



fo(co) 



VIA| 2 -(M ; 



I 2 
A 



+ t 



4T^| A p-(^)2 



(12) 
(13) 



3. Tunneling self-energy 

One dimensional chains with an average intra-chain 
hopping T are supposed to be connected by a weak in- 
terchain hopping t±. Interchain hopping corresponds to 
c t c a, encoding the destruction of a spinless fermion at 
site "a" and its creation at site a in a neighboring chain, 
as well as to the reversed process. The fermion creation 
operator is decomposed in the right (R)- and left (L)- 
moving components R and L according to 



1 

71 



"a,R 



-a,L 



with 



J a,R 



e ikFXa XR~{x a ) and c+ 



(14) 



Xl{Xa) 



(site "a" is at coordinate x a along the chain) . The fields 
Xr(l){ x ) are slowly varying as a function of the coordi- 
nate x along the chain. The identity 



leads to the self-energy 



1 1 
1 1 



(15) 



(16) 



for hopping from a to a, where the entries correspond to 
different left and right components. 



C. Weak localization-like loops 



2. Ballistic Green's functions 

Ballistic propagation within a given CDW chain is de- 
scribed by Eq. The CDW Green's functions take the 
form 



9 a >M 



A B 
C D 



with 



.4 



g (uj) exp {ik F {x a - xp)) 



(7) 
(8) 



A weak localization loop in a normal metal is shown on 
Fig. [5^. The generalization to subgap transport^ involves 
a Hikami bo x 59 ' 60 containing an "advanced-advanced" 
or a "retarded-retarded" transmission mode^i. Fig. [2}d 
shows a weak localization-like diffuson self-crossing in a 
3D gapped system (such as a superconductor^!) below 
the gap energy, and the same type of diagram from higher 
order terms in the tunnel amplitudes in normal metal- 
insulator-normal metal (NIN) structures are discussed in 
the forthcoming Sec. IHIl The generalization to quasi-lD 
CDWs (see Sec.HV[) corresponds to Fig. and a process 
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FIG. 3: (Color online.) Schematic representation of (a) a tight-binding normal metal-insulator-normal metal (NIN) tunnel 
junction with a hole pierced by a magnetic flux $ in the insulator (b) the fD NIN tunnel junction—, and (c) a geometry with 
an insulating ring pierced by a flux $, and (d) the resulting h/2e oscillations of the conductance with a negative magnetocon- 
ductance at low field. The notations a, a, /3 and b for labeling the interfaces are shown on the figure. 



of cotunneling through a CDW is shown on Fig. [2ji. As 
a direct consequence of Eq. (j4]), the transmission mode 
associated to Fig. [5Jd and Fig. [2t is /i/2e-periodic as a 
function of the flux $ enclosed in the loop (see the dis- 
cussion in Sees. IIIII and |IV[) . The magnetic field phase 
factor [see Eq. (TJJ] accumulated by the "advanced" and 
"retarded" Green's functions cancel with each other in 
the diagram on Fig. |2Ji, which is thus not modulated by 
a magnetic field. 



III. TUNNELING THROUGH A BAND 
INSULATOR 

We start with a tunnel junction consisting of a band 
insulator inserted in between two normal electrodes. The 
assumption on the dispersion relation (see Fig. [T] and 
Sec. Ill A[) is representative of extrema of the dispersion 
relation at wave- vectors p n , with 2ir/p n much smaller 
than the thickness of the insulating layer traversed by 
tunneling electrons. 



Ring geometry within microscopic Green's 
function 



1. Transport formula 

We start with the Caroli, Combescot, Nozieres and 
Saint- Jamesi expression of the conductance Q{uj) of a 
ID N a IN fc junction: 



(17) 



471 " 2 "^" Pa,a{^)ta,aGa,p{^)tl3,bPb,b{^)tb,pGp ta {u)t a ^ 

where fku is equal to the bias voltage energy eV; p a , a {u) 
and pbfi{u)) denote the density of states in electrodes "a" 
and "&" at energy ftw; t a = t a>a = t a>a and t b = = tpj, 
are the hopping amplitudes for the contact with the 
points a and (3 at the extremities of the insulator to 
their counterparts a and b in the normal electrodes (see 
Figs. [3^ for planar interfaces, and Fig. [HJd for the ID 



modeli). Eq. ([17]) is valid to all orders in the tunnel 
amplitudes t a and tb, from tunnel junctions to highly 
transparent interfaces. The fully dressed advanced and 
retarded Green's functions G A p(u>) and G^ a {u}) describe 
propagation from a to (3 and from (3 to a respectively 
while including all possible excursions in the normal elec- 
trodes, as opposed to the notations g A Ju) and gp a {uj) 
restricted to electrodes isolated from each other, with 
t a = tf, = 0. Physically, electron propagation across the 
insulator from electrode N a at time r a to electrode Nb 
at time t& [corresponding to the advanced Green's func- 
tion G^ o(lo)] occurs within the same process as prop- 
agation of an electron backward in time from electrode 
N& at time to electrode N a at time T a [correspond- 
ing to the retarded Green's function G§ a (u)}. As we 
show, both electrons propagating forward and backward 
in time can bounce independently at the interfaces dur- 
ing a tunneling process, corresponding to contributions 
to the conductance of higher order in tunnel amplitudes. 
This amounts to evaluating a transition probability as 
the square of a transition amplitude, which contains in- 
terference terms corresponding among others to weak 
localization-like subgap tunneling. 



2. Expansion of the transport formula 

Considering a geometry in which a quasi- ID ring made 
of a band insulator is connected to two normal electrodes 
(see Fig.|3t), the fully dressed advanced Green's function 
G^ 3 [see Eqs. (8) and (9) in Ref. Q 



G 



a,f3 



n A 

9a, /3 



]A (d' b ) A ' 



(18) 



with d A 



1 — 9a,a^O!,a9 a , a ta,a & n d d A 



9f,fi t fi,b9b,b t b,l3^ (d' a ) A - S^aWSoUVai ( d 'b) A = 

g A gt/3 : b9 A b tb,/3, is expanded in multiple crossings of the 



ring: 



G 



a,/3 



n A 
9a.fi 

d A d A 



(19) 
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FIG. 4: (Color online.) Notations used in Sec. IIIIBI for 
describing an insulating ring connected to two normal elec- 
trodes. 



The expansion in Eq. (|19p is justified by the damping of 
the electron wave-functions in the insulator^i: the sec- 
ond term with three crossings of the insulator in Eq. (fl9|) 
is small compared to the first term involving a single 
crossing if we suppose that half of the perimeter of the 
ring is large compared to the coherence length. Sep- 
arating propagation from a to (3 along the upper and 
lower branches of the ring and including the phase fac- 
tors related to the enclosed flux <£>, the combination of 
Eq. fT7 |) -([ig ]) to Eqs. JX^-JMl) leads to dominant h/2e 
oscillations of the conductance with the magnetic flux 
<fr, and to a negative magnetoresistance at low field (see 
Fig. EH). 

More precisely, the magnetic field is introduced by sub- 
stituting Eq. into Eq. (flU)) and keeping the lowest 
order oscillating term in the combination pGp a (uj) 
entering the expression of the conductance Giyi) given 
by Eq. (|17[) . Such term contains three advanced and 
one retarded Green's functions, each contributing to a 
factor exp (z7r<i>/20o). As a result of averaging over the 



Fermi phase factors (see Appendix |A")) . the lowest or- 
der oscillating part of the conductance is proportional to 
cos (47r$/(/)o) = Rc[cxp (4i7r<i>/0o)], with /i/2e-periodic 
oscillations of the conductance with the flux $, typical 
of a weak localization-like phenomenon. The sign of the 
cos (47r$/(/>o) oscillations in the conductance corresponds 
to a negative magnetoconductance at low field, and it 
is thus inverted as compared to standard weak localiza- 
tion in a normal metal (see Fig. [5] in the forthcoming 
Sec. WM- 

Similar features were obtained by Latyshev et aL— in 
experiments on CDWs (see the items 3. and 4. in the 
summary of experiments in the Introduction). The co- 
incidence is explained in Sec. IIVI by the same underly- 
ing weak localization-like subgap tunneling processes for 
band insulators and CDWs. 

B. Wave function approach 



We discuss now magneto-oscillations of the tunneling 
current via wave function matching 6 -^. The goal is to 
confirm qualitatively the Green's function approach in 
Sec. IIII A| and to present concrete results for the mag- 
netic oscillations. Let us consider the geometry on Fig. [4] 
in which a ring (made of an insulator in which wave 
functions are damped over a length £ and oscillate with 
wave- vector kp) is connected to two normal electrodes 
(with propagation of plane waves). The band insulator 
is supposed to have the same BCS dispersion relation 
as described above (see Sec. Ill A|) . The corresponding 
wave-vectors ±Uf ± l eac » to oscillations of the wave- 
function with period Xp = 2n/kp and to a damping over 
the insulator coherence length £. Following Refs. I62ll63l . 
the wave-functions iPl{x) in the left normal electrodes, 
ipiijji) in the upper branch of the ring (of length li), 
^2 (2/2) in the lower branch (of length I2) and iPr{x') in 
the right normal metal are given by 



^1(2/1) 



^2(2/2) 



cxp(-ikpx) + bexp(ikpx) 

. f, a , ; exp [-i^-yxAA sin (k F (h - y x )) exp ( Vl — — ) 
sin (Mi) V 4>o J \ ? / 

. tf , ; exp (i^-{h - yi)A ) sin(fc F ?/i)exp (-^ ] 
sin (Mi) \ <Po J V ? / 

. 7T~~r\ ex P I - i -[-y2Ao ) sin (k F (l 2 - y%)) exp ( — 
sm(M2) V 4>o J V 

tpgi ( 2ir . \ . /— 2/2 

■ n 1 \ ex P % ~T '2 -2/2) A) sin %2 sin —— 
sm (k F l 2 ) V 90 J \ t, 

b' e~xjp(ikpx'), 



(20) 



'y% - h 



(21) 



where we choose a radial gauge with a component Aq of correspond to the transmission and backscattering ampli- 
the potential vector, and where the coefficients b' and b 
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tudes, and ip a , ipp, ip a ' and ippi denote the amplitudes of 
the components of the insulating wave-function localized 
on the different nodes. Wave- function matching takes the 
form Vl(0) - Mh) = MO), Vte(O) = ^i(O) = Mh), 
and matching of the derivative of the wave-function at 
the left and right nodes is given by 



dipL 
dx 

dx 



(0) 
(0) 



dtp 



dip 2 



9 l W ■> 
ayi ay 2 

S'/'l , dtp; 



9y 



i(0) 



dy 2 



(0) 
(h) 



(22) 
(23) 



as deduced from integrating the Schrodinger equation 
with respect to coordinates over a small area including a 
node. The transmission coefficient 



CO 

o 



10 



~i 1 1 

R/^=0.5 ■ 
R/^=1 • 
R/^=2.5 ▼ 




T ($) = |6'($)| 2 (24) 

is then evaluated by averaging over the Fermi phase fac- 
tors UfIi and kph, and it is shown on Fig. [5]for different 
values of The matching equations (|20H23[) reduce to 
those of a normal metal 62,63 for R/£ -C 1 and a positive 
magnetoconductance with (po/2 periodicity is then recov- 
ered in this range of (as for standard weak localiza- 
tion in a normal metal where increasing a magnetic field 
suppresses localization). As seen from Fig. [SI the magne- 
toconductance is almost </>o/4-periodic at the cross-over 
R ~ £ and becomes again (/> /2-periodic for R/£ ^> 1, 
with a shape of oscillations characteristic of the evanes- 
cent wave-function weak localization-like oscillations dis- 
cussed above in Sec. IIII Al within microscopic Green's 
functions. We obtain /i/2e-periodicity for R/£ ^> 1, but 
not exactly of the form cos (Aircp / cpo) obtained above in 
Sec. IIII A 21 within microscopic Green's functions. This 
is because of a remaining positive magnetoconductance 
within a small low field region around 6 = rupa/2 (with 
n and integer), in agreement with Ref. |64|. Higher order 
harmonics play also a role because of the highly trans- 
parent contacts used for the ring geometry on Fig. 2] [see 
Eqs. ([22]) and |23|) ]. 



FIG. 5: (Color online.) Transmission coefficient of the ring 
made of a band insulator with the dispersion relation on Fig.fj] 
as a function of the flux <£> enclosed in the loop, normalized 
to the flux quantum <po = h/e. Different curves correspond 
to the different values of R/£ shown on the figure. The value 
of \f = 2n /kF is small compared to £. 



the conductance, and involve a transfer of 2kp momen- 
tum from one interface to the other across the insulating 
electrode. In a quasi-lD geometry, electrons above the 
gap on the right branch (label R) of the BCS-like disper- 
sion relation can propagate physically to the left or to the 
right according to their group velocity (see Fig. [TJ. This 
means that a recoil of the insulator is a necessary con- 
dition for weak localization-like subgap tunneling with 
quasi-ID insulators in the geometry of Fig. [6] Moreover, 
transfers of 2k f momentum are the hall-mark of CDW 
Andreev processe o 12 ! 13 ' 14 , which suggests a connection 
between the special type of tunnel junctions considered 
here and the CDW case (see Sec. IIVj). 

Finally, disorder in the normal electrode tends to lo- 
calize the electron wave functions in the vicinity of the 
interfaces. The processes on Fig. [5] are thus facilitated 
by diffusive motion in the normal electrodes. 



C. Coupling to a momentum channel 



Now, we note that propagation across the insulator 
defines a "tunnel" of cross section area £q, with £o ~ 
ooCf/A the coherence length (ao is the lattice spacing 
and ep the Fermi energy). Such a narrow channel is 
compatible with the following additional assumption: the 
insulator (with the dispersion relation on Fig. [T]) consists 
of linear ID chains perpendicular to the interfaces. Elec- 
trons on the right and left branches (denoted by R and V) 
are taken into account according to Sec. Ill Al We deduce 
that the specific set of R and L labels on Fig. [5^ does not 
contribute to the conductance once the summation over 
the Fermi oscillations in different channels is carried out. 
By contrast, branch crossing at the interfaces (see the R 
and L labels on Fig. [5Jd) contribute for a finite value to 



IV. SUBGAP TUNNELING BY A QUANTUM 

INTERFERENCE EFFECT AROUND A 
NANOHOLE IN A CHARGE DENSITY WAVE 

Now, we present our main result and show that weak 
localization-like subgap tunneling discussed above for 
tunnel junctions (Sec. IIII[) leads to a quantum interfer- 
ence effect also in charge density waves, therefore explain- 
ing the experiment by Latyshev et al.— on oscillations of 
the CDW current around a nanohole. Before discussing 
weak localization-like subgap tunneling in Sec. IIVBI we 
provide in Sec. IIV Al a mechanism of transport around a 
nanohole in a CDW. The disordered case is discussed in 
Appendix [B] 
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FIG. 6: (Color online.) Weak localization-like loops in a nor- 
mal metal-insulator-normal metal junction with two sets of 
right-left labels. The right-left branches and therefore the mo- 
mentum of the tunneling electrons is conserved in (a). The 
conductance is vanishingly small for (a) because of the ab- 
sence of propagation from a to /3 of the advanced-advanced 
transmission mode [g£ ^ R R g£ R R ~ - see Eq. @] . 
(b) shows a process leading to a finite tunneling current 
\da R n9a p l l i s limited by the insulator coherence length 
- see Eq. p)l]. Transfers of momentum by Ap — ±2kF prop- 
agate across the insulator according to the arrows, in parallel 
to evanescent wave charge tunneling. 



A. Coupling between the sliding motion and weak 
localization-like subgap tunneling 

1. Notations 

The mechanism discussed now for CDW transport, 
based on weak localization-like subgap tunneling around 
a nanohole, is first discussed for a few interrupted chains 
coupled by a transverse hopping t_i_. The more realistic 
case of a normal island around the nanohole will be dis- 
cussed afterwards. From the point of view of notations, 
the interrupted chains are labeled by fc — 2, ...,N — 1, 
and are connected by transverse hopping terms to the 
two uninterrupted chains k — 1 on top and k — N on 
bottom (see Fig. [7] for N = 3 and Fig. rjfc for N = 5). 



1 are interrupted at positions x 

(k) 



(fc) 



Chains k = '. 

at the left of the nanohole, and at xf 1 at the right of the 
nanohole (see x^ and x^ 1 on Figs.JTjD and [5] for a single 
interrupted chain). 



2. Deceleration and acceleration of the sliding motion in 
interrupted chains 

Assuming an overall sliding motion, the phase (fk(x, r) 
at position x along chains k = 1 and k — N and time r 



follows difk(x, t)/8t = ujq, with loq the sliding frequency, 
corresponding to the relation (see for instance Ref. I33T ) 



j (fc W) = - 



7r dr 



(25) 



between the CDW current and the time derivative of the 
CDW phase variable in chain k. The boundary condi- 
tion on the hole leads to d(pk( x T T ) I '9t = at x — x^ 
(i = 0,1 and k = 2,...,N— 1) (no collective current 

is flowing across x^ and x± in the direction parallel 
to the chains), and to dipk(x,T)/dr — ujq for x <C x^ 

(k) 

and x 3> x\ (the collective sliding motion is recov- 
ered far away from the nanohole). The resulting pro- 
file of d(fik(x,T)/dT (see Fig. [7J: for N = 3) corresponds 
to the conversion of the CDW current into a normal 
curren t 65 i 66 i 67 i 68 dpk(x,T)/dr, emitted according to the 
arrows on Fig. [7Jd, from the intermediate chain labeled 
by fc: 



dpk(xQ k) ,T) _ _e 5Vfc(xcTVr) 
dr 7r dxdr 

as deduced from the continuity equation. 

3. Charge accumulation 



(k) 



(26) 



Quasiparticles emitted from the slowing down of the 
sliding motion at the left of the nanohole are reabsorbed 
at its right where the sliding motion accelerates, leading 
to charge accumulation at the extremities of the inter- 
rupted chains at the left of the hole, described by the 
chemical potential Sp, > A. The current Ik(®) emitted 
from chain fc is given by 



-y 



T k ^ m (<S>,t ± ,hLu)d(hLj), (27) 



where 7fc^ m ($, t±, huj) is the total dimensionlcss trans- 
mission coefficient at energy fku transfering electrons 
from chain fc at the left of the nanohole to chain m at its 
right (see chains fc and m on Fig. ^ . The transmission 
coefficient Tk^ m in Eq. (j27|) is a multichannel general- 
ization of 7o in Eq. Similarly to the case of band 
insulators (see Sec. the subgap tunneling current for 
the processes on Fig. |S]is /i/2e-periodic as a function of 
the enclosed flux $ [see Eq. (|28[) below obtained from 
Eq. 0]. The oscillations appear only in the presence of 
the sliding motion (item 2. in the Introduction). The 
smallness of interchain couplings in the transmission co- 
efficient can be balanced by the integral over energy in 
Eq. (|2T|) . up to the large value of the Peierls gap in the 
compound NbSe3 used by Latyshev et al&. The non 
modulated part of the current for a single interrupted 
chain is proportional to Iq <~ (e//i)(tj_/T) 4 (<5/i — A) while 
the modulated part I m od = ( e /^)(^J-/7 1 ) 8 A is indepen- 
dent on Sfi-A. The ratio I mod /Io = (t//T) 4 A/(Sp- A) 
is thus a function of the values of t±/T and of (Sfi— A) /A. 



8 



nanohole 




chain (3) f*^^ 
/ 



chain (2) 

\ 
chain (1) 



FIG. 7: (Color online.) Schematic representation of (a) a film of NbSe3 pierced by nanoholes with a magnetic field B along the 
2 axis. The current is supposed to flow on average along the x axis parallel to the chains, (b) shows a nanohole interrupting a 
single CDW chain along the x axis. The arrows on (b) represent schematically the emission and absorption of normal carriers 
due to the slowing down and acceleration of the sliding motion at the left and right of the nanohole respectively. The profile 
of dtp n (x,T)/dr along chains n= f,2,3 is shown schematically on (c). 




chains (1) and (3) 
chain (2) 



chain (3) 
chain (2) 

chain (1) 




FIG. 8: (Color online.) The same weak localization-like tun- 
neling process as on Fig. [2J; for a nanohole in a CDW, as on 

Fig.m 



However, it was already noticed** that the damaged 
region of the CDW around a nanohole is likely to be 
normal, in which case the relative amplitude of the mod- 
ulation may be large, as suggested by Fig. [5] for a single 
channel. The corresponding tunneling process is shown 
on Fig. [5] at energy smaller than the normal region level 
spacing 5. Propagation across the normal region sup- 
ports "advanced-advanced" transmission modes because 
the level spacing 6 plays the role of a gap. The value 
of 6 is comparable to the Peierls gap, as it can be seen 
from the estimate Hvf/D, with D the diameter of the 
nanohole, leading to the rough estimate (5/fcg ~ 100 K, 
with vp ~ lO^iis -1 as an order of magnitude of the 
Fermi velocity. As expected, the level spacing of an ob- 
ject of size D ~ £ is comparable to the Peierls gap A. 



Nanohole 



Normal region 
around the nanohole 




Advanced-retarded 
transmission mode 

Advanced-advanced 
transmission mode 



R 

A 



FIG. 9: (Color online.) Schematic representation of a weak 
localization-like tunneling events involving an "advanced- 
advanced" transmission mode, and an "advanced-retarded" 
transmission mode in the presence of a normal region around 
the nanohole, at energies smaller than the island level spacing. 



B. Weak localization-like subgap tunneling 
transmission coefficient 



We evaluate now the weak localization-like subgap 
tunneling transmission coefficient given by the diagram 
on Fig. [51 in the absence of normal region around the 
nanohole as on Fig. [9l It takes the form 

T ($, t ± , fiw) = T (M H(fiw) cos , (28) 

The factor H(/kj) encodes the damping of subgap trans- 
mission: 



zz(huj) = exp 



(2) (2)' 



CM 

and T (Hu>) takes the form 



exp 



(2) (2) ' 

x l x 



CM 



(29) 



^(M - Re 4 | 3 oH| 4 {( 5 oM) 2 + (/oM) 2 } (30) 



{boMI 2 + l/oMI 2 )} 
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with g (uj) and f (u) given by Eqs. (Q2J) and (13]). The 
term \g Q (w) | 2 + |/o(w)| 2 corresponds to the "advanced- 
retarded" transmission mode in the lower branch (see 
Fig. [7H) and the term (g (uj)) 2 + (f (id)) 2 accounts for 
the "advanced-advanced" transmission mode in the up- 
per branch. The later involves transmission of momen- 
tum in parallel to evanescent wave tunneling, similarly 
to the tunnel junction discussed in Sec. [TTT] (see Fig. [BJd). 
The coherence length at energy Hlo is given by 

£(w) = hvF/\/\A\ 2 — (fuo) 2 , with vf the Fermi velocity. 

The sliding motion induces a time dephasing of 
weak localization-like tunneling because the CDW phase 
evolves in time in the course of weak localization wind- 
ing. Such couplings are however not probed in exper- 
iments such as in Ref. @ because the sliding motion is 
slow compared to the time scale H/A, with A the Peierls 
gap. 

A ballistic evaluation of the transmission coefficient 
in the CDW chains is justified by the fact that tunnel- 
ing quasiparticles travel over extremely short distances 
corresponding to the diameter of the hole, compara- 
ble to the CDW coherence length of order 10 nm in 
experiments^. Similar results were obtained in a differ- 
ent limit by treating disorder in the ladder approximation 
along Ref. [6l| (the principle of the calculation is detailed 
in Appendix |B|) . 

V. CONCLUSIONS 

To conclude, we discussed for CDWs a mechanism of 
tunneling via quantum interference effect initially pro- 
posed for superconducting hybrida^"^. The considered 
tunneling mechanism combined to charge accumulation 
due to the deceleration of the CDW at the approach 
of the nanohole leads to the same features as in the 
experiment by Latyshev et al. 6 on h/2e oscillations of 
the CDW current around a nanohole (see the Introduc- 
tion): 1. Weak localization-like loops due to higher order 
terms in the tunnel amplitudes are present even without 
nanoholcs. They induce no oscillations in the magnetore- 
sistance in the absence of nanohole because of the absence 
of charge accumulation in this case; 2. No charge accu- 
mulation is present with nanoholes but without sliding 
motion; 3. /i/2e-periodic oscillations of the resistance 
as a function of the magnetic flux are obtained with 
nanoholes and with sliding motion because the diameter 
of the nanohole is comparable to the coherence length, so 
that weak localization-like loops enclose approximately 
the same area as the nanohole; 4. The positive magne- 
toresistance at low field is already obtained for the nor- 
mal tunnel junction; 5. The only energy/temperature 
scale in weak localization-like subgap tunneling is the gap 
or the level spacing of the metallic island, comparable in 
magnitude to the Peierls gap. 

An underlying issue is whether the experiment by 
Latyshev et al£ provides evidence for an interference 
effect associated to the collective quantum mechanical 



CDW ground state. The collective momentum channel 
may be realized by a recoil in the specific case of a quasi- 
1D insulator. The analogy between CDWs and quasi-lD 
insulators and on the other hand the expected normal 
island in the Latyshev et al. experiment 5 shows that 
propagation through the CDW condensate is not a nec- 
essary condition for the modulations of the resistance as 
a function of a magnetic field. 
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APPENDIX A: TUNNELING THROUGH AN 
INSULATOR WITH PLANAR INTERFACES 

1. Multichannel transport formula 

The exact generalization of the ID transport formula 
(see Eq. (|17p) to a multichannel N a IN(, tunnel junction 
with extended interfaces (see Fig. [3^) is given by 

X Pb^M^.^GgaJw)^,^, (Al) 

where an underlying tight-binding lattice is assumed. 
The density of states p ai .aj connects the two sites dj and 
Oj at the same interface. The summations over and 
(k,l) run over all sites at the N a I and IN& interfaces re- 
spectively. The transport formula given by Eq. (|A1|) is 
represented on Fig. [TOb , where the bold red and dashed 
blue lines correspond to p and G^ a respectively. The 
densities of states in the normal electrodes are repre- 
sented by the connection of a wavy line on the figure. 

2. Perturbative expansion 

Higher order tunnel processes containing weak 
localization-like contributions are obtained by expanding 
systematically Eq. (|A1|) in the tunnel amplitudes accord- 
ing to the Dyson equations. Fig. [TO] b-f shows combina- 
tions of some terms in G^ p to some terms in Gp The 
weak localization-like diagrams on Fig. [TOB and e are 
made of three "advanced" and one "retarded" Green's 
functions, and look similar to the diagram on Fig. ^Bp 
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(d) (e) (f) 

FIG. 10: (Color online.) Top view of Fig. [3^. showing the first terms in perturbation theory in the tunnel amplitude for a 
normal metal-insulator-normal metal (NIN) tunnel junction. The exact transport formula is shown schematically on (a), where 
the thick solid (dashed) lines correspond to the fully dressed advanced (retarded) Green's functions. The wavy lines indicate 
the insertion of the density of states in the normal electrodes (see Fig. |3p). The other panels correspond to a few of the first 
low order terms in the perturbative expansion in the tunnel amplitudes. As for usual weak localization in a normal metal (see 
Fig. [2^,), the weak localization-like tunneling processes on (d) and (e) lead to dominant h/2e oscillations as a function of the 
magnetic flux once the summation over all channels has been carried out. 



corresponding to a disorder self-energy in the bulk of a 
superconductor. 



3. Averaging over the conduction channels 

The diagrams appearing in perturbation theory are av- 
eraged over the different "channels" in real space, corre- 
sponding to a summation over the discrete tight-binding 
sites at which diagrams cross the interfaces. The result of 
the channel averaging procedure depends on assumptions 
about the insulator band structure. 

Averaging over the conduction channels amounts to 
evaluating the integrals 



A 



A 



g A (R a ,f3,oj)g R (R a:l 3,uj) 
2^ 



Kp I A(„ . .\„R 



1 / Kp N 

2 \kFRa,/3 



(A2) 

g A (r,Lu)g H (r,Lu)dr (A3) 
{hw) 2 _ / 2R aJ} 



■ exp 



g A {R a , l3l uj)g A (R a ^,uj), 



|(A4) 
(A5) 



where g A (R a ^,oj) is given by Eq. ([T]), and where the 
last equality is valid for fku < A below the gap. To 
summarize, Aharonov-Bohm like oscillations are washed 
out by channel averaging with the band structure on 
Fig. [1] and /i/2e-periodic weak localization- like diagrams 
on Figs. [TOH and e contribute to leading order (in the 
tunnel amplitudes) to the oscillations of the conductance 
as a function of the magnetic flux. The same conclusion 
holds for the ring geometry in Sec. IIII A 21 



+ 



+ 



R 



R 


/ 


V 


R 






A 












I 


I 


+ 


R 











+ ... 



FIG. 11: (Color online.) The ladder series for the transmis- 
sion coefficient. The vertical blue lines correspond to impu- 
rities in the ladder approximation. The green boxes corre- 
sponding to Eq. (|B5fl indicate the recursive elimination of the 
impurities from one extremity of the laddepSi. 



APPENDIX B: EVALUATION OF THE 
TRANSMISSION COEFFICIENT OF A 
DISORDERED CDW 

1. Green's function of a disordered CDW 

The Green's functioni^ of a disordered CDW is eval- 
uated in the Born approximation as in the superconduct- 
ing case, where we introduce forward- and backward- 
scattering potentials u and v. With the notation in Ref. 
l6ll in the superconducting case, we find 



+ £,t 3 + Afi 
(fiu) 2 - f - A 2 ' 



(Bl) 



with fi and T3 the Pauli matrices given below, with ui = 
u(l + a), £ = £ + [5 and A = A(l + 7), with 



1 



1 



TfVA 2 -UJ 2 T£,\/A 2 — bJ 2 



(B2) 
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(' 



u\ 2 + \v\ 2 ) I x/ k 



1 



7 



u\ 2 + \v\ 2 ) I X F 



\u\ 2 fx/' 



2tt (fi^) 2 - a 2 - e 

dk £ 

2n (Huj) 2 - A 2 - £ 2 



(B3) 
(B4) 



2tt (So;) 2 - A 2 - £ 2 ' 



where Ai? is the Fermi wave- vector, ftw the energy, A 
the CDW gap and £ the kinetic energy with respect to 
the Fermi level. The notations 77 and t& stand for the 
forward and backward scattering times respectively. We 
keep in the following calculations a finite shift (3 of the 
chemical potential that does however not enter the prop- 
erties that we consider. 



2. Evaluation of ladder diagrams for the 
transmission coefficient 



To include disorder—, the transmission coefficients in 
the ladder approximation (see Fig. [TT|) are evaluated ac- 
cording to Smith and Ambegaokar- 1 by iterations of 



F(r n ) 



(B5) 



where u, v and w are Gaussian distributed random vari- 
able. The Pauli matrices f n are such that 











TO = 


i 







«T 2 

We find 



1 

-1 



, T 3 



1 

x 1 

1 

-1 



(B6) 
(B7) 



F(f ) = (~A + B)(l~a)hujT 3 (B8) 

-3,4(1 + y)Aif 2 
F(fi) = (A-B)(l+ 7 )Af 3 (B9) 

—ZA(1 — a)hujiT2 
F(ir 2 ) = 3(A + B)(l + 7 )Afo (BIO) 

—3A(1 — a)Hu)Ti 
F(t 3 ) = (-A-B)(l-a)tUJTQ (Bll) 

-A(l + 7)Afi, 



with 



A = 



(8k 2 F + 19q 2 )\u\ 2 



F (A 2 (l+ 7 ) 2 -/iu; 2 (l-a 



,2^1/ 



(Sfcg. + lQg 2 )!^! 2 



(B14) 



Acting twice with F according to G = F 2 leads to the 
closed 2x2 equations 



G 



'T2 
T3 



c d \ ( if 2 
d' c> { T 3 



(B15) 
(B16) 



with 



a = (A 2 -B 2 )(l -a) 2 (M 2 (B17) 

-9A(A + J B)(1 + 7 ) 2 A 2 
b = A(10A-B)(l-a)(l+j)HuA (B18) 
6' = -(A + B)(10A- B)(l -a)(l+7) (B19) 

x/iwA 

a' = 9^ 2 (1 - a) 2 (ftw) 2 (B20) 
- J 4(A-5)(1+ 7 ) 2 A 2 , 

and expressions of the same type for c, d, c', d': 

c = -3(A- B)(2A + B)(1 -a){l +7) (B21) 
xfiwA 

d = -9A[( J 4 + 5)(1 + 7 ) 2 A 2 (B22) 
-a) 2 (M 2 ] 

c' = (A — B) [(A + B)(l - a) 2 (hhj) 2 (B23) 

-A(1 + 7 ) 2 A 2 ] 
d' = ZA{2A + B){1 -a)(l + 7 )/iu;A. (B24) 

The final step is to decompose the initial condition on the 
eigenvectors of G and evaluate the coefficients in matrix 
geometric series such as 



+00 



E 6 " 

n=l 
+00 

E 6 " 



n=l 



0) -^(o)+*(ftl "-^ 



We find 



M = X^-^ + (B27) 
B 1 = X^-i^ + Y^-^ 2) (B28) 



1 — A-i 



A,, = -A'— K -±^^+X' 



B 2 = -A 



1 - A_ 

A_ 
1-A 



(B29) 



'j^-^+X'^-^, (B30) 



with 



x = I (1 - b '- a 

2\ {V - a) 2 + Aa'b 
2\ {V - a) 2 + Aa'b 



(B31) 
(B32) 
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■gi = (B33) ITl both sectors leads to an expression of all right-left 

(b' — a) 2 + Aa'b components of the transmission coefficient of a disor- 

1 / i- — \ dered CDW. The later can be used to evaluate the weak 

* = 2V a+ ± v( a - & ) + 4a &J (B34) localization-like subgap tunneling diagrams. 

, (-21 b' - b ± V(a - b') 2 + Aa'b 
^ = Vl 2& ' (B35) 

and ip^/ — 1. Carrying out the same calculation in the 
sector (172,73) and evaluating geometric series like 

+00 / » \ 

E 6 " # (o (B36) 
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